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Abstract. The cellular automata with local permutation invariance are 
considered. We show that in the two-state case the set of such automata 
coincides with the generalized Game of Life family. We count the number 
of equivalence classes of the rules under consideration with respect to 
permutations of states. This reduced number of rules can be efficiently 
generated in many practical cases by our C program. Since a cellular 
automaton is a combination of a local rule and a lattice, we consider also 
maximally symmetric two-dimensional lattices. In addition, we present 
the results of compatibility analysis of several rules from the Life family. 

1 Introduction 

The number of possible local rules for g-state cellular automaton defined on a 
(k + l)-cell neighborhood is double exponential of k, namely, q qk+1 . It is natural 
to restrict our attention to special classes of local rules. S. Wolfram showed \I\ 
that even simplest 2-state 3-cell automata, which he terms elementary cellular 
automata (there are 2 2 = 256 possible rules for these automata), demonstrate 
the unexpectedly complex behavior. 

We consider here a class of symmetric local rules defined on a (k + l)-cell 
neighborhood. Here by a 'symmetry' we mean a symmetry with respect to all 
permutations of k cells (points or vertices) surrounding (fc + l)th cell, which time 
evolution the local rule determines. The reasons to distinguish such rules are 

— The number of possible symmetric rules is the single exponential of polyno- 
mial of k (see formula (J2J below). For example, for k = 8, q — 2 (the case 
of the Conway's game of Life) the number of symmetric rules is 262144 « 
2.6 x 10 5 , whereas the number of all possible rules sa 1.3 x 10 154 . 

— The symmetry of the neighborhood under permutations is in a certain sense 
a discrete analog of general local diffeomorphism invariance which is be- 
lieved must hold for any fundamental physical theory based on continuum 
spacetime. 1 

1 The symmetric group of any finite or infinite set M is often denoted by Sym(M). If 
M is a manifold, then a diffeomorphism of M is nothing but a special — continuous 
and different iable — permutation from Sym(M). 



— This class of rules contains such widely known automata as Conway 's Life. 
In fact, as we show below, any symmetric rule is a natural generalization of 
the Life rule. 



2 Symmetric Local Rules and Generalized Life 
2.1 Symmetric Rules 

We interpret a (k + l)-cell neighborhood of a cellular automaton as a k-star 
graph, i.e., rooted tree of height 1 with k leaves. We call this the k-valent neigh- 
borhood. We adopt the convention that the leaves are indexed by the numbers 
1, 2, . . . , k and the root is numbered by k + 1. For example, the trivalent neigh- 
borhood looks like this 



A local rule is a function specifying one time step evolution of the state of root 



We consider the set Rs k of local rules symmetric with respect to the group 
of all permutations of leaves, i.e., variables x\, . . . , x k - We will consider also the 
subset Rs k+1 C Rs k of rules symmetric with respect to permutations of all k+ 1 
points of the neigborhood. For brevity we shall use the terms k-symmetry and 
(k + 1)- symmetry, respectively. 

Obviously the total numbers of k- and (fc + l)-symmetric rules are, respec- 
tively, 



2.2 Life Family 

The "Life family" is a set of 2-dimensional, binary cellular automata similar to 
Conway's Life [2], which rule is defined on 9-cell (3x3) Moore neighborhood 
and is described as follows. A cell is "born" if it has exactly 3 alive neighbors, 
"survives" if it has 2 or 3 such neighbors, and dies otherwise. This rule is sym- 
bolized in terms of the "birth" / "survival" lists as B3/S23. Another examples of 
automata from this family are HighLife (the rule B36/S23), and Day & Night 
(the rule B3678/S34678). The site |3] contains collection of more than twenty 
rules from the Life family with Java applet to run these rules and descriptions 
of their behavior. 
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Generalizing this type of local rules, we define a k-valent Life rule as a 
binary rule on a fc-valent neighborhood, described by two arbitrary subsets of 
the set {0, 1, . . . , k} . These subsets B, S C {0,1, ... ,k} contain conditions for 
the Xk — > x' k transitions of the forms — > 1 and 1 — > 1, respectively. Since the 
number of subsets of any finite set A is 2> A ', the number of rules defined by 
two sets B and S is equal to 2 fc+1 x 2 k+1 = 2 2k+2 , which in turn is equal to 
|2"| evaluated at q = 2. On the other hand, different pairs B/S define different 
rules. 

Thus, we have the obvious 
Proposition. For any k the set of k-symmetric binary rules coincides with the 
set of k-valent Life rules. 

This proposition implies, in particular, that one can always express any sym- 
metric binary rule in terms of "birth" / "survival" lists. 



2.3 Equivalence with Respect to Permutations of States 

Exploiting the symmetry with respect to renaming of q states of cellular au- 
tomata allows us to reduce the number of rules to consider. Namely, it suffices 
to consider only orbits (equivalence classes) of the rules under q\ permutations 
forming the group S q . For counting orbits of a finite group G acting on a set R 
{rules, in our context) there is the formula called Burnside's lemma. This lemma 
states (see., e.g., @]) that the number of orbits, denoted \R/G\, is equal to the 
average number of points R 9 C R fixed by elements g G G: 

Thus, the problem is reduced to finding the sets of fixed points. 

Since we are mainly interested here in the binary automata, we shall consider 
further the case q = 2 only. For this case the combinatorics is rather simple. 

Specializing and © for q — 2 we have the numbers of binary k- and 
(fc + 1)— symmetric rules, respectively, 

N Sk = 2 2k+2 , (5) 
^Vs fc+1 = 2 fe + 2 . (6) 

The group S2 contains two elements e = (0)(1) and c = (01) (in cyclic 
notation). After S. Wolfram, we shall call the permutation c G S2 "black-white" 
(shortly BW) transformation. 

Since the number of fixed points for e is either (JSJ or JBJ , all we need is to 
count the number of fixed points for the BW transformation in both k- and (k + 
l)-symmetry cases. The rules in both these cases can be represented, respectively, 
by the bit strings 

aia 2 ■ ■ • a 2 k+2 (7) 

and 

aia 2 ---a k+2 (8) 

in accordance with the tables 
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One can see from these tables that the BW transformation acts similarly on 
both strings Q and JSJl, namely, 



«1«2 ' ' -OLn-XO-r, 



BW 



a n a. n _ x ■ ■ ■ a 2 a 1 , 



where bar means complementary transformation of bit, i.e., a = a + 1 mod 2, 
and n = 2k + 2 or n = k + 2. 
The fixed point condition 



aia 2 ■ ■ ■ a n -\u r , 



a n a n -\ ■ ■ ■ a 2 a\ 



(9) 



implies that the string is defined by a half of bits, i.e., by k + 1 or by (k + 2)/2 
bits. In other words, the numbers of different bit strings satisfying condition 10 
are, respectively, 2 k+1 and 2^ k+2 ^ 2 (for k even). For the (k + l)-symmetry case 
with odd k — 2m + 1 condition @ leads to the contradiction 



a 



m+l 



o 



m+l j 



which means zero number of bit strings in this case. 

Summarizing our calculations (recollecting that G = S 2 in formula |0J, i.e., 
\G\ = 2), we have: 



— in the case of fc-symmetry 

• number of BW-symmetric rules 

Ns k BW 

• number of non-equivalent rules 

N S k /BW = 

— in the case of (k + l)-symmetry 

• number of BW-symmetric rules 



■ife+i 



,2fc+l 



+ 2\ 



2 k/2+l i£k= 

if k = 2m - 



1. 



• number of non-equivalent rules 

f 2 fe+1 + 2 fe / 2 if k = 2m, 
1 2 fc+1 if k = 2m - 



Ns k+1 /BW 



(10) 
(11) 

(12) 
(13) 



For example, for trivalent rules the numbers are 

Ns 3 bw = 16, Ns 3 /bw = 136, N Ss+1 bw = 0, 



Ns 3+1 /bw = 16. 
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3 Assembling Neighborhoods into Regular Lattices 

fc-valent neighborhoods can be gathered into a lattice in many ways. In fact, any 
fc-regular graph may serve as a space for a cellular automaton. In applications 
a cellular automaton acts, as a rule, on a lattice embedded in a metric space, 
usually 2- or 3-dimensional Euclidean space. Any graph, being 1-dimensional 
simplicial complex, can be embedded into E 3 . When dealing with symmetric local 
rules it is natural to consider equidistant regular systems of cells. We consider 
here only two-dimensional case as more simple and suitable for visualization of 
the automaton behavior. Most symmetric 2D lattices correspond to the tilings 
by congruent regular polygons. The number of different types of such tilings is 
rather restricted (see, e.g., 0). To denote a fc-valent lattice composed of regular 
p-gons we use the Schlafli symbol {p, k}. 

3.1 2D Euclidean Metric 

There are only three regular lattices in E 2 (see Fig.^l. 



3-valent, {6, 3} ^-valent, {4, 4} 6-valent, {3, 6} 




Fig. 1. All regular lattices in E 2 



Since real computers have finite memory, cellular automata can be simulated 
only on a finite lattice. Usually the universe of a cellular automaton is a rectangle 
instead of an infinite plane. There are different ways to handle the edges of the 
rectangle. One possible method is to fix states of the border cells. This breaks 
the symmetry of the lattice and is thus not interesting for us. Another way is to 
glue together the opposite edges of the rectangle. 

The toroidal arrangement is a standard practice, but it would be interesting 
to study cellular automata on nonorientable surfaces also. There are 3 different 
identifications of opposite sides of a rectangle: the torus T 2 , the Klein bottle K 2 , 
and the projective plane P 2 . All these spaces (the gluing does not affect their 
Euclidean metric) are shown in the figure below. 
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We need to check whether the regular lattices like in Fig.^ i- e -j with the Schlafli 
symbols {p, k} = {6, 3}, {4, 4}, {3, 6} can be embedded in T 2 , K 2 , or P 2 . To do 
this we must solve the system of equations 

V-E + F = x(M), P F = kV = 2E, (14) 

where V, E, F are numbers of vertices, edges, and faces, respectively, %(M) is the 
Eulcr characteristic of a manifold M. Since x(T 2 ) = x(^ 2 ) = an d x(P 2 ) = 1, 
we see that the regular lattices are possible only in the torus and the Klein bottle 
and impossible in the projective plane, as well as in any other closed surfaces. 

Summarizing, for Euclidean metric there are only 3-valent hexagonal, 4- 
valent square and 6-valent triangular regular lattices in E 2 , T 2 , and K 2 . 



3.2 Hyperbolic Plane H 2 

The hyperbolic (Lobachevsky) plane H 2 allows infinitely many regular lattices. 
Poincare proved that regular tilings {p, k} of H 2 exist for any p, k > 3 satisfying 

p ' k ^ 2 - 

For example, let us consider the octivalent Moore neighborhood used in the 
Life family and shown in the figure below. 



The Moore neighborhood can not form regular lattice in the Euclidean plane 
since the distances of surrounding cells from the center are different. But in the 
hyperbolic plane regular 8-valent lattices exist and there are infinitely many of 
them. The simplest one is shown (using the Poincare disc model projection) in 
Fig. 12 




Fig. 2. Octivalent regular lattice {3,8} in H 2 
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3.3 Sphere S 2 



All regular lattices in the two-dimensional sphere correspond to the Platonic 
solids which are shown in Fig. [3J 

3-valent 

, A « 

{3,3} - tetrahedron {4,3} - hexahedron {5,3} - dodecahedron 




Fig. 3. All regular lattices in S 2 



In the sphere there are infinitely many other lattices, which are close to 
regular and may serve as spaces for 3-valent symmetric automata. They are 
called fullerenes. 

The fullerenes were first discovered in carbon chemistry in 1985, and this 
discovery was rewarded with the 1996 Nobel Prize in Chemistry. A model of 
the first revealed fullcrcnc — the carbon molecule Ceo — is displayed in Fig. 
0] (the figure is borrowed from Later there were discovered other forms of 
large carbon molecules with structural properties of fullerenes (larger spherical 
fullerenes, carbon nanotubes, graphenes). Their unique properties promise they 
will have an important role in future technology in particular, in nanotechnology 
engineering. 

From a mathematical viewpoint, the structure of fullerene is a 3-valent convex 
polyhedron with pentagonal and hexagonal faces. In terms of graphs, fullerene 
can be defined as a 3-regular (3-valent) planar, or equivalently, embeddable in 
M = § 2 , graph with all faces of size 5 or 6. 

Let us generalize slightly this definition assuming that M is not necessarily 
§ 2 , but may be closed surface of other type, orientable or nonorientable. Then 
the Euler-Poincare equation together with the edge balance relations gives the 
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Fig. 4. A model of Ceo carbon molecule (buckyball). 



system of equations 

V-E + f 5 + f 6 = X (M), W = 5/ 5 + 6/ 6 = 2E, (15) 

where f$ and fe are numbers of pentagons and hexagons, respectively. The gen- 
eral solution of this system is 

h = 6x(M), (16) 
V = 2/ a + 10 X (M), (17) 
E = 3/ 6 + 15 X (M). (18) 

We see that generalized fullerenes are possible only in the sphere § 2 , in the pro- 
jective plane P 2 , and in the torus T 2 and Klein bottle K 2 . Attempts to consider 
surfaces with greater genus lead due to (|16fl to senseless negative numbers of 
pentagons. Thus, we have for all generalized fullerenes: 

^ = 2/ 6 + 20, E = 3/ 6 + 30, / 5 = 12, sphere § 2 ; 

V = 2f 6 + 10, E = 3/ 6 + 15, f 5 = 6, projective plane P 2 ; 

V = 2/6, E = 3/6, h = 0, torus T 2 , Klein bottle K 2 . 

We see that any fullcrcnc in § 2 or in P 2 contains exactly 12 or 6 pentagons, 
respectively, and arbitrary number of hexagons. In the case of torus or Klein 
bottle the fullerene structure degenerates into purely hexagonal lattice without 
pentagons considered already in subsection 13.11 Note that in carbon chemistry 
one-layer graphite sheets (and similar structures with few pentagons or hep- 
tagons added) are called graphenes. 

4 Canonical Decomposition of Some Rules from Life 
Family 

In this section we present the canonical decompositions of relations for three rules 
from the Life family: the standard Conway's Life, HighLife and Day & Night. 

Recall that canonical decomposition |7] of relation on a set of points is the 
representation of the relation as a combination of its projections onto subsets 



of points. This decomposition is discrete analog of compatibility analysis of al- 
gebraic and differential equations. To apply our approach we interpret the local 
rule x'g — f (xi,X2, ■ ■ ■ , Xg) as a relation on 10 points x±, X2, ■ ■ ■ , xg, x' Q . 

The HighLife is interesting since for it the replicator — a self-reproducing 
pattern — is known explicitly. For Conway 's Life, the existence of replicators 
is proved, but no example is known. 

The name Day&Night reflects the BW-symmetry of this rule. Perhaps, it 
is the first (conceivably found "by hand" ) automaton from the Life family with 
this symmetry property. Note that there are exactly 512 such automata in the 
Life family, and they all are generated easily (for time < 1 sec) by the C program 
mentioned in [7j. 

It might be easier to grasp the relations, if they are written in the form of 
polynomials over the field F2. In the below formulas we use elementary symmetric 
polynomials in variables x±,...,xs- Recall that the kth (of degree k) elementary 
symmetric polynomial of n variables x\, . . . , x n is defined by the formula: 

0~k (xi , ■ ■ ■ , X n ) ^ Xi 1 Xi 2 ' • • Xi k . 

I<ii<i2<---<ik<n 

Hereafter we use the notations Ok = Ok ■ ■ ■ , %s), a\ = <7fc (x%, . . . , Xi, . . . , xs) , 
u\ 3 = Ufc (x%, . . . , x~i, . . . , £j, . . . , xs) ■ The indices i, j, k, I satisfy l<i<j<fc< 
I < 8. 

The polynomials representing Conway's Life, HighLife and Day&Night 

take the forms, respectively 

-Pconway's Life = A + Xg {<7 7 + (T 6 + <7 3 + (7 2 } + (Jj + <J 3 , (19) 
-PHighLife = x'g + Xg {cr 3 + <7 2 } + ^6 + CT 3 , (20) 
PDay&Night = Xg + Xg {<7 7 + (J 6 + (T 5 + er 4 } + (T 8 + <7 7 + (J 6 + (T 3 . (21) 

Note that these polynomials have degrees 8, 6, 8 and numbers of terms 185, 169, 
256, respectively, so application of the Grobner basis technique to their analysis 
may take some time (about 1 hour on 1.8GHz AMD Athlon notebook with the 
Maple 9 Grobner procedure). Our program computes the decompositions for 
time < 1 sec. 

The decompositions are: 

— Conway's Life: 

x'g {0-3 + 0-2 + 1} + 0-7 + 0-3 = 0, (22) 
x'gxg {a\ + a\} + x'g {a' 2 + l} + x 9 {a\ + a\ + a\ + a\) = 0, (23) 

4{4 + 4 + ^ + l} = 0, (24) 

x'g (xg + 1) {of + a l 2 J + of + l} = 0, (25) 

XgXiXjXkXi — 0. (26) 
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HighLife 

Xg {03 + <7 2 + 1} + (7 7 + 03 = 0, (27) 

x' 9 x 9 {a\ + a\) + x' 9 {4 + a\ + l} + x 9 {4 + a\ + 4 + 4} = 0, (28) 

x' 9 {<t\ + <r\ + 4 + a\ + 1} = 0, (29) 

+4{4 i +^+^'+4'+4' + ^' + l}=0, (30) 

XgXg^XjXfcX; = 0. (31) 

— Day&Night 

Xq {(77 + CT 6 + (75 + 0"4 + 1} + cr 8 + (T 7 + (T 6 + (T 3 = 0, (32) 

ZgZg {4 + 4 + a\ + 4} +4 {4 + 4 + 4 + °i + *} 

+:r 9 {4 + 4}+4=0, (33) 

4 {<# • 4' • 4' • 4 j • 4' • i} • < o. (34) 

Note that system (I34[l of prime relations can be combined into the reducible 
relation (for terminology like prime and reducible see \7\) 

4 {4 + 4 + 1} + 4 = 0, 

which looks nicer (the polynomial representation of relations is somewhat 
artificial), but depends on larger set of variables. 

We see that the decomposition for Day&Night differs essentially from the de- 
compositions for Conway's Life and HighLife having resembling structures. 

5 Conclusions 

We proved that the cellular automata from the Life family are nothing but binary 
automata with the local rules symmetric with respect to all permutations of the 
outer cells in the neighborhood. 

Then we showed that the number of non-equivalent with respect to renaming 
of states fc-valent symmetric binary local rules is equal to 2 2k+1 + 2 k . Considering 
the 3-valent case — the next step up after the 2-valent Wolfram's elementary 
automata — we see that the total number of non-equivalent symmetric rules is 
only 136. 

All interesting 3-valent 2-dimensional regular (or almost regular) lattices are: 

— hexagonal lattice {6,3} in the plane E 2 , in the torus T 2 and in the Klein 
bottle K 2 

— tetrahedron {3, 3}, hexahedron (cube) {4, 3} and dodecahedron {5, 3} in the 
sphere § 2 
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— fulleren.es in S 2 and in the projective plane P 2 . 

Combining these lattices with 136 symmetric rules we obtain a class of cellular 
automata which looks like quite available for systematic study. Since these au- 
tomata have more interesting geometry than the elementary ones we may expect 
more interesting behavior of them. 
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